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Abstract 

The methods of non-equihbrium thermodynamics of systems with an interface 
have been applied to the study of thermionic emission processes in abrupt 
semiconductor junctions, including the effects of surface states . Our analysis 
covers different situations of interest, concerning unipolar and bipolar systems, 
either degenerate or not, with or without interface states and under steady and 
non-steady conditions. In this way, a complete phenomenological modelling of 
thermionic emission-interface states processes has been proposed, overcoming 
some of the inherent limititations of the currently existing models. 
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I. INTRODUCTION 



Metal-semiconductor contacts have been extensively studied in the literature |T|], 0. At 
the beginning, basically two different approaches were followed to describe their rectifying 
properties: the diffusion approach [Q, in which the main current limiting mechanism was the 
diffusion of carriers through the semiconductor; and the thermionic emission (TE) approach 

in which the dominant process was the thermal emission of carriers over the contact bar- 
rier. These two approaches were synthesized later |Q, giving rise to the so-called thermionic 
emission- diffusion (TE-D) models, in which both current limiting processes were taken into 
account. A similar development occurred in the more recent area of heterojunctions, in 
which both TE [§, and TE-D models were also proposed. 

These models have been further generalised by including the effects of interface states 
that very often are present at semiconductor junctions, see for instance, P], 0, [|T0| . 
One may call these generalizations thermionic emission-interface state-diffusion (TE-IS-D) 
models. In some cases, mostly in the context of metal-semiconductor (MS) junctions, a thin 
insulating layer was introduced to account for some non-ideal properties [||], [Q. Its main 
effect was to allow the electric potential to be discontinuous across the junction. Therefore, 
in this paper we will distinguish between abrupt junctions (continuous electric potential) 
and non-abrupt ones (discontinuous electric potential). For many purposes, the models 
mentioned above suffice to describe most of the properties of semiconductor junctions 

Despite their wide application, a general phenomenological derivation of those models 
is still lacking in the current literature. This may be due to the fact that the TE-IS and 
diffusion processes are, a priori, of a different nature, which suggests applying to the former a 
microscopic-kinetic approach, while to the latter a phenomenological one. As a consequence, 
this fact has caused a situation where the models are not formulated in a general enough 
way, in the sense that some situations or processes have not been explored, mainly those 
with respect to the TE-IS processes. 
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It is our purpose in this paper to present a complete analysis of the TEl-IS processes in 
abrupt semiconductor junctions. This study overcomes some of the inherent limitations of 
the previous models and provides a proper framework to study these junctions in a more 
complete way. 

We would like to draw special attention to the fact that we will develop the whole anal- 
ysis in the same phenomenological framework as where diffusion processes may be treated, 



111 , allowing us then to present a unified derivation of the TE-IS-D models in abrupt semi- 



conductor junctions. In this way all the situations and processes may be studied, including 
the effects of degenerancy, bipolarity, non-stationarity or interface state filling, which are 
treated only in an approximate way in most of the models mentioned above. 

The organisation of the paper is as follows: In section II, we summarize the main results 
concerning non-equilibrium thermodynamic description of semiconductor junctions. In this 
framework, and with the help of the formalism of the internal degrees of freedom, we derive 
in section III the TE law for abrupt unipolar junctions without interface states. This result 
is further generalised in section IV, to include the effects of interface states or bipolarity. 
Finally, in section V we summarize our main results. 



II. NON-EQUILIBRIUM THERMODYNAMIC DESCRIPTION OF 
SEMICONDUCTOR JUNCTIONS 

When two different semiconductors (or, in general, two different materials) are in con- 
tact, a very narrow region exists in between, which exhibits peculiar macroscopic behaviour. 
It has been shown |[T^, that the properties of this region can be explained in a simple 
way, by substituting it by a two-dimensional interface, with its own macroscopic properties, 
and by assuming that the bulks to will behave as in the absence of the junction just up to this 
interface. This introduces a discontinuous description of the system, in which the relevant 
fields are the extrapoled bulk fields (discontinuous through the interface) and the surface 
fields (to describe surface properties), |TT| JT^. The discontinuity in the bulk fields makes it 
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necessary to introduce some discontinuity relationships in the description of these systems, 
and the presence of the two-dimensional interface the proper surface equations to describe 
its properties. A systematic and thermodynamically consistent way to derive a complete 
set of equations (bulk, interface and discontinuity equations) to describe these systems, is 



provided by the theory of non-equilibrium thermodynamics of systems with interfaces [1TB 
|T8| . Precisely, by following this approach, a one dimensional phenomenological model for 
drift-diffusion processes in abrupt semiconductor junctions was derived in which in- 



cluded, apart from the usual drift-diffusion bulk equations, the corresponding surface and 
discontinuity relations. 

Of particular interest for the present paper are the discontinuity relationships for the 
drift-diffusion processes, for they describe the TE-IS processes that take place at the semi- 



conductor junction. To arrive at these discontinuity relations, we first evaluated in |TT| the 
surface entropy production corresponding to a semiconductor junction when only drift and 
diffusion processes are considered. For a one-dimensional system one obtains ||11||, 



a' = -^J^-X^-^J--X~ (1) 
for the case in which interface states are present, and 

= -\;J^ ■ [$\^ (2) 

when they are not. In Eqs. ([|) and (H), we have defined = (Jj^, ^2"); where stands 
for the perpendicular component to the interface of the current density number of electrons 
(k=l) and holes (k=2), evaluated at the left (— ) and right (+) hand sides of the junction. 
Similarly, we have X± = (Xf , Xf) , = [^J-), with X+ = X^ = 

and [(pk]- = 4>k ~ 4>k ^ with 0^'* being the corresponding electrochemical potentials (quasi- 
Fermi levels) for the k-th carrier in the bulks (±), or at the interface (s), defined as 

<Pt = f^t + QkV^, 0^ = ^± + g4y]+; k=l,2 (3) 

where fj,^'^ are the chemical potentials (Fermi levels), qk the charge of the k-th carrier, and 
y± the electric potential with [V]+ = i(\/+ + l^"). All the magnitudes defined above are 



evaluated at the interface, that is, at x = 0. For abrupt junctions, the electric potential is 
continuous across the junction, then one has [V]^ = 0. 



/^From these entropy productions the following discontinuty relations may be derived [jTT 

1 



J^ = --F^{X\X~) 



(4) 



which is valid when interface states are present. Otherwise, the corresponding relation is 

J-" = -^m]-) (5) 

In the above equations, and F are some unspecified functions satisfying the equilibrium 
condition 



F±(0,0) = F(0) = 

the second principle of thermodynamics 

X+ ■ F+ + X- ■ F- >0 
• F > 

and the Onsager's symmetry relations 



(6) 



(7) 
(8) 



dF^ 




dF^ 




k,J = l,2 




eq 


dX^ 


a, /9 = +, - 

eq 


dFi 










eq d[<P,]^ 


eq 



(9) 



(10) 

Finally, it should be noted that the functions F^ (F) may also depend on the bulk number 
densities, (surface number density n|) and on temperature, T ITTl . 

For the case under consideration, in which recombination-generation or tunnelling pro- 
cesses are neglected, one usually assumes that the transport of current across the interface 
is controlled by the thermal excitation of carriers over the barrier (or barriers) that appears 
at the junction. For this reason, we will call Eqs. (^) and (|^) the TE-IS and TE relations, 
respectively and the search for explicit expressions for them is the main purpose of next 
sections. 



III. THERMIONIC EMISSION RELATION FOR UNIPOLAR ABRUPT 

JUNCTIONS 

As commented previously, the derivation of all these expressions will be carried out by 
means of the methods of non-equilibrium thermodynamics. Due to the intrinsic non-linear 



nature of the emission processes, one has to go beyond the usual linear phenomenological 
relations, and propose more general ones. One way to systematically derive such a non-linear 
relations, still under the framework of non-equilibrium thermodynamics, is provided by the 
formalism of the internal degrees of freedom. This will be the approach we will follow in this 
paper. For the sake of clarity, we will devote the next subsection to introduce the internal 
degrees of freedom in the context of semiconductor junctions, postponing the derivation of 
the TE relation for next subsection. 

A. Internal degrees of freedom and semiconductor junctions 

In the thermodynamic description of a system one has to distinguish between transport 
and relaxation processes. For the former, transport equations accounting for the evolution 
of quantities, as density, velocity, etc. can be formulated starting from the conservation 
laws of mass, momentum, etc. Transport equations then describe the evolution of those 
quantities, essentialy due to two mechanisms namely, diffusion and convection. Relaxation 
processes involve transformations ocurring at each point of the system whose mechanism 
is intrinsically non-linear and is usually modelled by the pass through a potential barrier. 
In non-equilibrium thermodynamics language this process may be viewed as diffusion in 
an internal space, the space of the internal variable. These scheme has been satisfactorily 
applied to a number of physical situations ||l3l, 



Concerning semiconductor junctions the following example will illustrate the way of 
introducing internal variables to arrive at a complete characterization of the system. Let 
us consider the emission of a carrier through the interface. When the carrier reaches the 
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interface, it has several options: to jump directly into the other side of the junction, or into 
an empty surface state, or simply to go back. The net effect of these processes is known to 



be non-linear [|rT|, in the sense that linear expressions for Eqs. (^) or (^), do not correctly 
describe them. To derive these non-linear relations within the framework of non-equilibrium 
thermodynamics, one may use the internal variables. In this alternative picture, one may 
consider that the different processes a carrier may go through when it reaches the interface, 
are well represented by a continuous diffusion process in the internal space, over an energy 
barrier which matches the discontinuity points at the interface. The carriers following this 
continuous process will be called activated carriers, and the coordinate that describes it , 
the internal coordinate. 

By introducing an internal space for each elementary process, whose number may be 
predicted directly from the expression of the entropy production, the net rate corresponding 
to them may be evaluated, allowing one to arrive at an explicit expression for Eqs. 
and d). 

It is worth emphasizing that both pictures (internal degrees of freedom and macroscopic 
discontinuity) leads to the same entropy production |]T2|], [|T^, and it is this fact that enables 
us to use the former to derive explicit non-linear constitutive relations, provided the latter 
one is unable to. 



B. Thermionic emission relation 

For the case of unipolar, i.e. those with only a single carrier, abrupt semiconductor 
junctions without interface states, the surface entropy production may be written as Eq. (H), 

^: = -^JM- (11) 

For the sake of definitess, we have taken the electrons as the single carrier and the cor- 
responding magnitudes have been denoted by the subindex e. This entropy production 
corresponds to a single elementary process, as can be seen by comparing it with the one 
corresponding to an unimolecular reaction of the form 

7 



provided one identifies with the reaction rate and [^J- with the affinity (see [|T^, [|T^). 
Furthermore, this equivalence between chemical reactions and emission processes, allows 
us to directly write the corresponding constitutive relation for the latter case by simply 
identifying the corresponding magnitudes in the former case. In this way we obtain, in the 
general case of a degenerate (non-ideal) system, [|l^, ||18| , 



J+ = ksl (l - e^[^='-) (13) 

with 

(14) 



where ks is Boltzman's constant, (3 = {ksT)^^, 7 is the internal coordinate, having values 
in the interval [71,72] of the internal space, D{'-)) is the diffusion coefficient in the internal 
space, f/(7) is the internal energy barrier, 72(7) is the density number of activated electrons at 
point 7, and /(7) = /(n(7), T) accounts for the degenerancy effects, and is defined through 
the state equation 

0(7) = ^BTln/(7)n(7) + t/(7) (15) 

where (/>(7) is the generalised chemical potential of the activated electrons. For the case of 
a non-degenerate system, / is only a function of T, which considerably simplifies Eqs. ([M]) 
and (pISl). In Eqs. (|13D and (|1^ we have used (^(71) = 0~ and 0(72) = (pt 1 fo^ these 



points the internal space makes contact with the external one. Moreover, in what follows, 
we will also identify /(n(7i),T) = f-{n-,T), /(n(72),T) = f+{n+,T) and f/(7i) = [/", 
U{j2) = U^- In the non-degenerate case f^'"''^ = (N^)^^, where A^^ is the effective density 
of states, which is only a function of T. As we have assumed the mass of the carriers to be 
the same on each side of the junction, we will then have that N~ = = Ne. 

A simplified expression for Eq. (|I1|) may be obtained by introducing a mobility into the 
internal space, /i, through Einstein's relation 
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= — d^) = ■ ^'^^ 

We then obtain 

I = T - (17) 

For the purposes of this paper, it is convenient to define a transition coefficient as follows 

A-+ = kB^^ . (18) 



With this definition, and by substituting Eq. ([T7|) in it, we have 

= ^/7^^/i(7)-ie/^(^W~^(^i))/(7)c/7 ■ ^^^^ 

Note that for an abrupt junction, f/(7) is a given function independent on the density 
number. Moreover fi may also be considered independent of the density number, even for 
a degenerate system. Hence, the main dependence of A_+ on the density is given through 
/. As a consequence, for a non-degenerate system, in which / is only a function of T, this 
transition coefficient will be independent on the densities. 

Similarly, we may define the inverse transition coefficient, A-|__, as 

A+_ = A_+e^^^ (20) 

where AU = U{^2) — U{ji) = — U~ represents the discontinuity in the electron energy, 
which for an abrupt junction is independent of the density. 
By using this coefficient we may rewrite Eq. (pTsD as 

J+ = X+.e-'f< [l - e^^'^^^-) (21) 



where e(f)^ = ~ 4>e ■ K is easily seen that Eq. pl|), or alternatively Eq. (|13|), satisfy 
the unipolar version of the thermodynamic requirements presented in section II, Eqs. 
and d^) ( no symmetry relation is needed here). Moreover, as predicted in that case, the 
discontinuity relation only depends on the thermodynamic force, in this case [4>e\-i and on 
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n^, through A+_, and T. On the other hand, for the non-degenerate case, A+_, is only a 
function of T, and we denote it by A^*^. 



Let us show that Eq. (^3|), in the non-degenerate case, is precisely equivalent to the 
boundary conditions used in the TE-D theories . By simply noting that for that case one 
has 



n- = N^e"^'^^-^^^ (22) 



we can rewrite Eq. (|2l|) as 

eJ: = e^{nl-nt) (23) 

where n+ = Nee~'^^'^b . If we identify vr = we then see that Eq. (p3D is precisely the 
boundary condition used in in the context of MS junctions. A kinetic calculation for 
it gives [Q f R = ^^]^, with A* being Anderson's modified constant. Therefore, we obtain 
is ,as mentioned above, only a function of T. 
Moreover, by using Eq. ( p2D again, we may rewrite Eq. (^) as 

eJt = {n:e-^^^ - n+) (24) 

which corresponds to the boundary condition proposed in |Q in the context of heterojunc- 
tions. It is worth pointing out that Eqs. (123|) and (0) reduce, respectively, to Bethe's 0] 
and Anderson's expressions, in the limit in which diffusion processes are neglected. 

IV. THERMIONIC EMISSION-INTERFACE STATE RELATIONS FOR ABRUPT 

SEMICONDUCTOR JUNCTIONS 

Our aim in this section is to generalise the results of the previous section, by including 
the effects of interfacial states or by considering bipolar systems, i.e. systems with two types 
of carriers. 
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A. Unipolar systems with interfacial states 



For this case the surface entropy production can be written as, see Eq. (0), 

- €) - - €) ■ (25) 

This expression is reminiscent of the entropy production corresponding to the set of in- 
dependent elementary processes concerning the interchange of electrons among bulks and 
surface, 



(26) 



labelled by a = 1, 2, 3, with reaction rates Ji = —J^ , J2 = — Je ^ ^3 = -Je ^ conjugated 
affinities Ai = X", A2 = X+, A3 = X~ + X+ = [0e]_. 



The equivalence with the chemical reaction case analyzed in [^], allows us to directly 

formulate the constitutive relations 

3 

Ji = Yl kB^iJa (1 - e"^-) z = 1, 2, 3 (27) 

a=l 

where z/j^ is the stoichiometric coefficient of the z-th electron ( we identify e~ = ei, = 62 
and e"*" = 63 ) participating in the a-th elementary reaction, A^ are the affinities defined 
above and 

g/3<)!'a(7c«i) 

^ ksPC, DM-'eP'^-^'-^U^) (1 + fi^) d7« ^^^^ 
where 0q(7qJ corresponds to the generalised chemical potential of the reactant'm the a — th 
reaction. As before, by introducing a mobility, jla-, through Einstein's relation we can rewrite 
Eq. ( p8|) as 

^ ^ p,:r^a{ia)-'eP^-^^'^)u^)d^^ ■ ^^^^ 
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For the elementary processes under consideration, Eq. pq), we have 



By substituting these values into Eq. 



1 1 

0-11 ■ (30) 
1 0-1^ 

we obtain the following constitutive relations, 



J7 = kf 



h l-e 



j; = k 



B 



+ h{l 

h (l-e^^^-)+/3 (l 



(31) 
(32) 



Note that the right hand sides of the previous equations are actually functions of the ther- 
modynamic forces, X^, [4>e]- (only two of them are independent because [4>e]- = X^ + X~ 
), and of T and nl. To see this last dependence, we note that although the phenomenological 
coefficients, at first sight, may also depend on the bulk densities, nf, this dependence gives 
rise to a dependence on and X^, through the equations of state 



(33) 



Furthermore, by construction, Eqs. (^Tj) and (^) satisfy the thermodynamic restrictions 
given through Eqs. (|]), (|^) and (H). We therfore conclude that they constitute acceptable 
thermodynamic constitutive relations. 

As in section III, a more useful form of these constitutive relations is obtained by defining 
the following transition coefficients 



A, 



I < J 
i > j 



(34) 
(35) 



where Zi = fiiii are the activities and Af/j,- = Ui — Uj. We then obtain a master-like equation 



i=3 



i=3 



i = 1,2,3 ; ij^j. 



(36) 



i=i i=i 

Note that for an abrupt junction AUij are given quantities, which allow us to write Eq. ( P3[ ) as 
a detailed balance -like expression, XijZ^'^ = XjiZ^'^ , where we have used e^^^'^ = e^^^^i 
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It should be pointed out that in the degenerate case, as can be seen from Eqs. (0) and (p5D, 
the transition coefficient may depend on the density number, and hence on the apphed bias. 
On the other hand, in the non-degenerate case, fa becomes independent on the density, 
and then the transition coefficients are only a function of T. In this case, we introduce new 
transition coefficients A°- = A^!^/""^, where as before we assume Z"'' to be independent on i 
and equal to Ng. Then the constitutive relations transform into |jl2| 



i=3 j=3 

J^ = ll>^%^3-ll>^l^^ ^ = 1,2,3 ; l^J (37) 
i=i i=i 

with X^juij^ = A^jTij^ and A°- only varying with T. 

Keeping these definitions in mind, we may rewrite Eqs. (|3TD and (^) as 

J+ = N,Xl,e-^^^^-*^^ (l - e^l^^l-) + iVeA°2e-^(^^^-^^) (l - e^^^^) (38) 
J- = NAlie-^^''^-'^^^ (l - e^^^") + NAlie-^^^'^-'f'^^ (l - e'^t'^^l") (39) 

which constitutes the generalization of Eq. (^1]), or 

J+ = XI, {n:e-^^^^-^^^ - nt) + A^, (^e'^^^-^^") - n+) (40) 
J; = A^i (n:e-^(^-^^") - nl) + A°, (n^e-^^^^"^^") - nt) (41) 



which also generalizes Eq. (|2^) . The corresponding expressions valid for the degenerate case 
may be obtained by simply substituting X^jN^. by Xij and by Zi. This simple rule then 
justifies the fact that in what follows only the non-degenerate expressions will be treated 
explicitly. 

Some final remarks are in order. First of all, it is worth noting that Eqs. (|38D and (p9D, 
or alternatively, Eqs. (|40|) and (|4l|), hold under non-steady conditions for which J+ ^ J~. 
Moreover, they also include all the possible transitions among interface states and bulks and 
all the situations of interest (degenerate and non-degenerate systems). Finally, it is worth 
noting that the number of parameters to be specified is relatively small: the transition 
coefficients A21, X'^,, A325 ^^e energy discontinuities U^' — U^, U^' — (clearly U~ — f/| 
is not an independent quantity). This makes the modelling of these processes not as difficult 
as it would seem a priori. 
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B. Steady-state conditions 



Under steady state conditions Eqs. (|38D - (^) can be simplified considerably. In this 



case one has that = , see and we then obtain 

J+ - M (\0 I I ^ e-^iU+-^e ) (l _ (AO) 

1 + |i e/5(t/e+-f/|) 
1 + ^e/3('^^-c/e^)e-/3[0e]- 



^32 



where, as noted in |11], the relevent variable under steady state conditions is [(pe]-- Alter- 



natively, we can also write 



The corresponding expressions for the degenerate case are obtained by using the transfor- 
mation rule formulated previously. 

It is worth mentioning that Eqs. P^ ) and (0) have a TE-like form, see Eqs (^TJ) and (P^, 
but they differ in the fact that now the prefactor incorporates the effects of the interfacial 
states. Furthermore, we have obtained the equations governing the population of those 
states under non-equilibrium conditions, Eqs. (^31) or (H), in their general form. From 
these relations, approximate expressions can be proposed. For instance, let us consider 
the case, very common in the current literature, in which it is assumed that the interface 
states are in permanent equilibrium with one of the bulk systems, say system — . This 
situation corresponds to the case in which the transitions e"^ — > e"*" are practically inhibited 
as compared to the — > e" transitions. In our formulation this means that 

A| = ^^-mi-ut) » 1 . (46) 

^23 A32 



Using this approximation in Eq. and considering moderate values for [0e]-7 we obtain 

€ ^ <p: (47) 
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which is the condition of permanent equihbrium between the interface and the — system, 
and which ultimately justifies the use of the bulk chemical potentials to control the filling 
of the interface states. Moreover, in this case, the prefactor in Eq. p^ ) reduces to Ag^ + A32, 
which is thought to introduce a small modification with respect to the commonly used TE- 
like expressions. It should be noted that for high values of [0e]- this approximation may 
not be justified, and non-trivial effects coming from the presence of interface states could 
appear. 



C. Bipolar systems 

In this subsection, we will consider the case of bipolar systems. For simplicity's sake we 
will start by considering the case in which interface states are not present. In this case the 
surface entropy production reads, Eq. (Q), 

< = -^[M-Jt-^[4>e]-Jt (48) 

where the subindex h stands for holes and e for electrons. It is easily seen that the elementary 
processes corresponding to this situation are simply 

e~ ^ e"*" 

h- . (49) 

Note that these processes are decoupled and hence they can be treated separately, proceeding 
along the lines indicated in section II. The constitutive relations for this case will then read 

J+ = iV,Af_e-^(^^^-<^=") (1 - e'^t'^^l-) = (50) 

Af_(n;e-'3(^^^-^^")-<) (51) 

J+ = NhX'^.e-^'^''^-'^'^^ (1 - e^l^'^l-) = (52) 

Xf_{n-,e-^^^^'^^^-nt) (53) 

where the meaning of the different magnitudes clearly follows from the notation. As before, 
the degenerate case is obtained by applying the corresponding rule. Again by construction 
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the corresponding thermodynamic restrictions are satisfied. Note that in this case, we have to 
specify four parameters, the transition coefficients X'^_, X^_, and the energy discontinuities 
f/+ — , Ujl — Uf^. For bipolar heteroj unctions the two carriers are hardly ever considered 
both together, for usually one of the discontinuities is much larger than the other inhibiting 
the crossing of the corresponding type of carrier [^. This fact is obviously recovered from 
our relations. 

The introduction of interfacial states in these systems does not complicate the situation 
very much, because due to the decoupling that exists between the two types of carriers, 
they can be treated separately, giving rise to a set of equations like the ones obtained in the 
previous subsection for each of the carriers. 



V. CONCLUSIONS 

In this paper we have presented a fully thermodynamic and general treatment of 
thermionic emission processes in semiconductor junctions. All the situations of interest, 
namely, for unipolar and bipolar systems, which can be degenerate or not, with or without 
interface states and under steady and non-steady conditions, have been analysed. 

The study has been carried out within the framework of non-equilibrium thermody- 
namics of systems with interfaces, incorporating the formalism of the internal degrees of 
freedom in order to obtain explicit expressions for the TE-I relations. This procedure has 
the advantage of allowing us a complete treatment of semiconductor junctions (bulks, inter- 
face and discontinuity relations) in a common framework which ensures by construction the 
thermodynamic consistency and completeness of the relations derived. 

Of special interest has been the derivation of general TE-IS relations to describe non- 
steady situations in the presence of interface states. From these expressions we have then 
derived the corresponding relations valid at the steady state, which have been shown to 
incorporate, in the general case, non-trivial contributions that come from the interface states, 
which might not be negligible in some circumstances. These interface state contributions 



16 



are under investigation for some junctions of interest. 

Finally, the implementation of the TE-IS expressions in the drift-diffusion model pre- 
sented in will provide one of the most general TE-IS-D models that has been proposed 
for abrupt semiconductor junctions, proving therefore the power of the application of the 
non-equilibrium thermodynamic methods to semiconductor systems. 
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